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The  report  presents  the  calculation  procedure  and  numerical  results  for  the  far  field  radiation  and  diffrac¬ 
tion  wave  patterns  around  a  ship  at  zero  forward  speed  The  calculation  procedure  uses  the  strip-theory  Ship 
Motions  Program  to  obtain  the  motion  transfer  coefficients  and  singularity  strengths  in  the  form  of  the 
integrated  Kochin  functions.  These  results  are  then  combined  with  an  asymptotic  expression  for  the  far  field 
wave  elevation  due  to  a  three-dimensional  oscillating  source  to  obtain  the  ship  wave  patterns 

The  wave  patterns  are  calculated  for  a  series  of  six  ships,  including  a  strut-like  hull,  excited  by  ambient 
waves  representative  of  Sea  States  2,  4,  6,  and  8.  The  radiation  waves  are  calculated  for  all  six  modes  of  ship 
oscillation.  The  diffraction  waves  are  calculated  for  wave  headings  ranging  from  following  seas  to  head  seas. 
The  accuracy  of  a  commonly  used  long  wavelength  approximation  to  the  Kochin  function  in  which  the  singular¬ 
ity  distribution  over  the  hull  surface  is  collapsed  to  a  line  distribution,  is  also  discussed,  y  , .  ri-  ■ 
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CALCULATION  OF  FAR  FIELD  RADIATION  AND 
DIFFRACTION  WAVE  PATTERNS  USING  THE  KOCHIN  FUNCTION 
APPROACH  FOR  A  SERIES  OF  SHIP  HULLS 


1.  INTRODUCTION 

The  calculation  of  the  force  and  moiton  coefficients  of  a  surface  ship  excited  by  ambient  ocean 
waves  is  a  problem  of  widespread  interest.  The  motion  coefficients  are  of  interest  to  ensure  the  safety 
and  comfort  of  the  crew  and  passengers,  as  well  as  the  stability  of  the  onboard  cargo.  The  force  coeffi¬ 
cients  are  of  interest  since  they  are  components  of  the  equations  of  motion,  and  also  enter  in  the  calcu¬ 
lations  of  the  structural  behavior  of  the  ship.  Practical  computational  approaches  may  be  broadly 
classed  in  two  categories:  strip  theory  where  interactions  between  transverse  ship  secuons  are 
neglected,  and  fully  three-dimensional  approaches.  An  example  of  the  strip  theory  approach  is  the 
DTNSRDC  Ship  Motion  Program,  SMP  [1.2]  while  an  example  of  the  second  approach  is  that  by 
Chang  [3],  A  comparison  of  these  two  approaches  for  a  series  of  ship  hulls  has  been  conducted  by 
Wang  and  Chang  (4j. 

A  related  proolem  of  considerable  interest  is  the  calculation  of  the  wave  profile  next  to  the  ship. 
This  is  of  interest,  for  example,  in  deck  wetness  and  slamming  calculations  Examples  of  such  calcula¬ 
tions  are  those  by  Lee  [51  and  Beck  [61. 

All  of  the  above  are  concerned  with  the  near  field  flow  and  wive  pattern  next  to  the  ship  hull. 
While  there  have  been  some  studies  outlining  the  general  calculation  procedure  for  the  far  field  wave 
pattern  [7,8],  actual  calculations  of  these  waves  for  typical  ship  hulls  appear  to  be  lacking.  This  report 
presents  the  calculation  procedure  and  numerical  results  for  the  far  field  radiation  and  diffraction  wave 
pattern  for  a  series  of  ship  hulls  at  zero  forward  speed  due  to  typical  ocean  waves.  Radiation  waves 
refer  to  those  caused  by  the  three  translational  and  three  rotational  oscillatory  motions  of  the  ship, 
while  diffraction  waves  are  due  to  the  scattering  of  the  incident  exciting  wave  by  the  ship  hull,  t2ken  to 
be  stationary. 

The  calculation  procedure  for  the  near  field  flow,  outlined  in  [8],  is  briefly  summarized  here.  The 
calculation  of  the  far  field  wave  patterns,  using  the  Kochin  function  approach,  is  discussed  in  greater 
detail.  It  is  shown  that  the  suggested  collapse  of  the  hull  surface  singularities  onto  a  line  [7,8]  is 
equivalent  to  an  approximate  calculation  of  the  Kochin  function,  assuming  i’nai  the  length  of  the  excit¬ 
ing  wave  is  much  larger  than  me  beam  and  draft  of  the  ship. 

Radiation  and  diffraction  wave  patterns  are  calculated  for  a  series  of  six  ship  hulls  ranging  in 
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length  from  140  to  990  ft,  four  wavelengths  of  the  exciting  waves  ranging  from  35  to  998  ft,  and  five 
wave  headings  ranging  from  following  to  head  seas.  The  wave  characteristics  are  given  in  terms  of  fig¬ 
ures  showing  the  variation  with  spatial  direction  of  the  amplitude  and  phase  ax  a  given  distance  from  a 
fixed  point  on  the  ship.  The  figures  show  typical  variations  of  the  wave  patterns  with  characteristics  of 
exciting  wave  and  hull.  The  accuracy  of  the  calculated  waves  using  the  approximate  Kochin  function 
procedure,  for  different  wavelengths,  is  also  shown.  The  report  concludes  with  a  summary  of  the  prin¬ 
cipal  findings. 

2.  DESCRIPTION  OF  CALCULATION  APPROACH 

Adopting  the  usual  assumption  of  Unearned  potential  flow,  the  problem  essentially  consists  of  two 
parrs:  the  evaluation  of  the  potential  flow  on  and  near  the  ship  hull  and  the  extrapolation  of  this  solu¬ 
tion  to  obtain  the  far  field  waves.  A  detailed  description  of  the  procedure  for  calculating  the  near  field 
flow  is  given  in  (81  and  hence  is  only  summarized  here.  The  evaluation  of  the  Kochin  function  and  its 
use  in  obtaining  the  far  field  wave  pattern  is  discussed  in  greater  detail. 

2.1  Calculation  of  Near  Field  Flow 

The  near  field  problem  consists  of  finding  a  potential  function  <£  which  satisfies  Laplace's  equa¬ 
tion,  the  free  surface  condition,  the  radiation  condition  of  outward  progressive  waves  at  infinity,  and 
the  kinematic  condition  of  no  flow  through  the  ship  hull.  In  the  commonly  used  Green’s  function 
approach,  the  procedure  is  to  place  a  series  of  singularities  on  the  hull  surface,  which  satisfy  the  first 
three  conditions,  and  determine  their  strengths  from  the  kinematic  condition  on  the  hull  surface. 
Largely  due  to  the  presence  of  the  free  surface  condition,  the  general  solution  of  this  problem,  even 
with  the  linearization  assumption,  is  sail  quite  complex. 

A  considerable  simplification  results  if  the  beam  and  draft  of  the  ship,  as  .'ell  as  length  of  the 
exciting  wave  are  all  taken  to  be  small  compared  to  the  length  of  the  ship.  In  th.s  case,  longitudinal 
interactions  between  ship  cross  sections  may  be  neglected,  and  the  complex  three-dimensional  problem 
reduces  to  a  series  of  simpler  two-dimensional  problems  for  various  transverse  sections,  or  strips.  An 
example  of  this  strip  theory  approach  is  the  Ship  Motion  Program  (SMP)  [1,2]  developed  at 
DTNSRDC.  For  the  zero  speed  case  considered  in  the  present  work,  it  is  shown  in  [4]  that  SMP  calcu¬ 
lates  force  coefficients  which  are  in  fair  agreement  with  those  gi’cu  by  *  three-dimensional  approach 
[3]  and  motion  coefficients  which  are  in  reasonable  agreement. 

The  coordinate  system  used  in  the  calculations  in  SMP  (and  also  used  in  the  present  work)  is 
shown  in  Figure  1.  The  origin  lies  on  the  undisturbed  free  surface,  at  the  longitudinal  center  of  grav¬ 
ity,  and  in  the  vertical  plane  of  symmetry.  The  x-aiis  is  directed  forward,  the  y-axis  to  port,  and  the 
z-axis  upwards.  In  this  coordinate  system,  following  and  head  seas  correspond  respectively  to  headings 
of  0  and  180  degrees.  It  should  be  noted  that  the  above  coordinate  system  and  sea  heading  convention 
do  not  correspond  to  those  used  for  input  and  output  in  SMP. 
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The  total  unsteady  potential  <f>  is  witten  as  the  following  sum 

4>  -  la  (fa  +  fa)  4-  £  e,m‘  (1) 

where  a  is  the  amplitude  of  the  incident  wave 

fa  ■  >  the  potential  of  the  incident  wave  of  unit  amplitude 

fa  is  the  difTraction  potential 

ij  is  the  amplitude  of  the  jth  oscillation 

4>j  is  the  radiation  potential  due  to  the  jth  oscillation  of  unit  amplitude 
co  is  the  circular  frequency. 

As  shown  in  Figure  l,  j  •»  1,3,5  correspond  to  the  motions  of  surge,  heave,  and  pitch  in  the  vertical 
plane,  and  j  -  2,4,6  correspond  to  the  lateral  motions  of  sway,  roll,  and  yaw. 

The  radiation  potentials  <t>j  are  obtained  by  placing  two-dimensional  Havelock  sources  on  a  series 
of  transverse  cross  sections  along  the  ship  longitudinal  axis.  The  strengths  of  these  singularities  are 
determined  by  the  kinematic  condition  of  no  flow  through  the  hull  cross  section  for  each  of  the  oscilla¬ 
tion  modes,  using  a  procedure  originally  implemented  by  Frank  [9],  It  is  shown  in  [7,81  that  the  dif¬ 
fraction  potential  fa,  which  is  obtained  from  the  kinematic  condition  of  no  flow  of  the  incident  wave 
through  the  stationary  ship  hull,  may  be  expressed  in  terms  of  the  heave  and  sway  radiation  potentials, 
as  follows 

fa  -  -  (.fa  -  isin  0  fa)e~‘laane  (2) 

where  0  is  the  wave  heading 

k  —  —  -  is  the  wave  number 
g  k 

to  is  the  circular  frequency  of  the  exciting  wave 
k  is  the  wavelength  of  the  exciting  wave. 

2.2  Calculation  of  Far  Field  Wave  Patterns 

To  fust  order,  the  strengths  of  the  sources  obtained  for  the  near  field  may  also  be  used  to  obtain 
the  far  field  waves  [71.  However,  the  sources  themselves  must  be  changed  from  those  for  two  dimen¬ 
sions,  Gjo,  to  those  for  three  dimensions,  Gjd-  For  example,  a  two-dimensional  source  gives  far  field 
waves  of  constant  amplitude  whereas  in  three  dimensions  they  decay  d\  t  to  cylindrical  spreading. 


Three  approaches  axe  outlined  in  [8)  for  using  the  near  field  singularity  strengths  to  obtain  the  fax  field 
waves:  by  using  the  damping  coefficients,  by  compressing  the  hull  surface  source  distribution  to  a  line 
distribution,  or  by  using  the  Kochin  Function  which  represents  an  integral  over  the  huii  surface  of  the 
source  strengths  weighted  by  an  exponential  function.  In  the  present  work,  the  Kochin  function 
approach  is  taken.  It  is  shown  below  that  the  line  distribution  approach  essentially  represents  a  long 
wavelength  approximation  of  the  Kochin  function. 

In  their  classic  and  comprehensive  survey  of  free  surface  flows,  Wehausen  and  Laitone  (10}  give 
the  following  asymptotic  expression  for  the  far  field  wave  elevation  £,  of  f.'e  radialion  wave  due  to  an 
oscillation  ^  of  unit  amplitude 

T-—H.  Gfae'-'  -  2 e'w*-w/4>  e™  (3) 

4sr  g  (  V  irR 

where  j  »  1,. .  .,6  denote  the  six  modes  of  ship  oscillation 

R  —  x2  +  y1  denotes  the  radius  in  the  horizontal  plane,  measured  from  the  origin  shown 
in  Figure  1 

0  —  tan-1  y/x  is  the  direction  in  the  horizontal  plane,  measured  from  the  x-axis 

Hj  is  the  complex  conjugate  of  the  Kochin  function  Hj 

G$d  is  the  asymptotic  expression  for  the  three-dimensional  source  potential  at  large  R. 

Several  equivalent  forms  for  Ht  are  given  in  (10],  In  the  present  work,  it  is  of  particular  interest  to 
consider  Hj  as  the  following  weighted  integral  of  the  singularity  strength  distribution  over  the  hull  sur¬ 
face  S 

Hj  U,  9) - jy  yj  (*,7 1,0  ("****•»&  (4) 

where  y  s  is  the  hull  surface  distribution  of  singularity  strengths  calculated  by  SMP  for  the  jth  mode  of 
ship  oscillation.  The  exponential  function  precisely  accounts  for  the  difference  in  the  wave  due  to  a 
source  located  away  from  the  origin.  Thus,  the  factor  ekl  gives  the  well  known  vertical  exponential 
damping  of  the  wavemaking  capability  of  a  source  situated  below  the  free  surface.  The  factor 
giUi  cos»+t)  «n  s)  gjves  change  in  phase  of  the  wave  originating  at  x  «  £,  y  —  i)  from  that  originating 
at  x  —  y  —  0. 

It  is  of  interest  to  note  that  the  far  field  asymptotic  expression  is  valid  at  relatively  short  distances 
from  a  given  source.  In  the  two-dimensional  case,  Marnyanskii  111]  shows  that  at  a  distance  of  0.3  A 


from  the  source,  there  is  Httle  effect  of  the  near  field  terms.  In  the  three-dimensional  case,  at  a  dis¬ 
tance  of  kj 6  from  an  oscillating  source  submerged  below  the  free  surface,  Liu  [12}  obtains  a  difference 
of  approximately  15%  in  the  wave  elevatioas  calculated  by  the  far  field  and  complete  formulas. 


For  the  actual  case  where  the  ship  oscillation  is  not  of  unit  amplitude  but  is  the  response  to  a 
wave  of  amplitude  a,  the  resultant  wave  elevation  Crj  must  be  written  in  terms  of  a  and  the  transfer 
function  7}  which  is  a  function  of  wave  frequency  w  and  wave  heading  /3 


{ kj  (*,  «.  R,fi)  -  Re 


aTj(w%0)  j  Hj  ( k,9 )  *'<**+•<-*/«>] 


~  Re  6,  R,  p)  ei<ut*-‘-*/4)) 


(5) 


^  /  k  _ 

where  Aj  (k,  9 ,  R,  fi)  —  a 7}  — -  /  — — —  Hj  (k,  9)  is  the  amplitude  coefficient  of  It  may  be 

£  V  iff 

noted  that,  for  the  sake  of  simplicity,  the  incoming  wave  is  here  assumed  to  have  zero  phase  with 
respect  to  the  origin. 


The  formula  for  the  resultant  diffraction  wave  £#7  is  conveniently  obtained  from  Eqs.  (2)  and  (5). 
Noting  that  the  ship  is  taken  to  be  stationary  and  hence  7}  does  not  enter,  £*7  is  given  by 

i/n  *  {k iv  +  imi 

(Hiy(k,  9 ,  p)  -  H7L  (*,«,  j8))  e,(W?+“'~ir/4,| 


-  Re  [ 


[A1V  (k,  8,  R,  p)  +  A 7t  (*,  9 ,  R,  p)e 


tikji  +»l— w/4> 


(6) 


where 


(vertical  plane)  and  antisymmetric  (lateral  plane)  components  of  Aly 
from  Eq.  (4)  with  y7r  —  — y3  and  y7t  “  /  sin  0  y2  e-*t 


are  respectively  the  symmetric 
and  Hly  and  H1L  are  evaluated 


2.3  Long  Wavelength  Approximation  of  the  Kochin  Function 

It  is  of  interest  to  investigate  the  consequence  of  evaluating  the  Kochin  function  by  assuming  the 
wavelength  X  to  be  long.  In  particular,  if  \  is  large  compared  to  the  transverse  dimensions  of  beam  and 
draft  of  the  ship,  the  exponential  terms  in  Eq.  (4)  may  be  approximated  to  first  order  by 

ek{  =  1  +  k{  (7a) 

e**’***"  :  „  ^  jin  ff,  1  sm  (if  yj  sin  y) 

=  l  +  iki)  sin  9.  (7b) 
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Considering  that  (£,  ij,  0  is  even  in  7;  for  the  vertical  modes  ij  —  1,3,5)  and  odd  m  n  tor  the 
lateral  modes  (j— 2,4,6),  the  integral  (4)  for  a  given  transverse  section  C  at  £  «*  £0  takc^  on  the  follow¬ 
ing’ two  fOnTiS  _ 

-  Jc>>  (^o.  h-  £)  (1  +  *0  (1)  dl  i-  1,3,5  (8a) 

Hj  =  -/At  sin  5  fcyj  (fo>y,0  (1  +  Ae£)  (tj)  dl  j  —  2,4,6  (8b) 

If  one  neglects  the  wave  clamping  terra  1  +  Eqs.  (8a)  and  (8b)  reduce  to 

Hj  =-Qj  J  -  1,3,5  (9a) 

Hj  =  — rArsin  9  M,  j  —  2,4,6  (9b) 


where  £>y  and  A f,  are  precisely  the  resultant  point  source  and  dipole  strengths  for  a  given  transverse 
section  C  given  in  181.  The  factor  ik  sm  9  appearing  in  Eq.  (9b)  is  accounted  for  in  [81  by  considering 


the  dipole  M )  to  create  waves  given  by  the  dipole  potential 


a  Gjp 

By 


(3) 


Using  the  definitions  given  in  Eq. 


dy^^R'^0'  U0) 

from  which  it  follows  that  at  large  values  of  R 

Gfo  -  ik  sin  9  GaiD  (11) 

Thus,  the  suggested  approach  of  collapsing  the  section  singularity  distribution  to  a  single  point  is 
equivalent  to  evaluating  the  Kochin  function  for  large  wavelengths,  neglecting  in  addition  the  vertical 
damping  tern  1  +  k(.  In  the  present  work,  this  terra  has  been  retained.  It  should  be  noted,  however, 
that  while  this  term  may  improve  the  approximation  at  large  values  of  X  (small  values  of  |A:{I),  it  may 
actually  worsen  it  at  intermediate  and  Large  values  of  I  Ac  4 ! .  For  example,  for  large  \k0,  l  +  Ar{  may 
become  negative  while  the  actual  minimum  value  of  ekl  is  0. 

3.  DESCRIPTION  OF  COMPUTER  RUNS 

3.1  Choice  of  Ship  Halls 

A  total  of  six  ship  hulls  were  considered.  The  principal  geometric  characterises  of  these  hulls  are 
shown  in  Table  1.  The  DE1006  destroyer  and  the  CVA-59  carrier  were  previously  investigated  in  the 
force  and  motion  coefficient  study  [4).  Their  lengths  are  respectively  308  and  990  ft.  In  order  to 
investigate  the  effect  of  appendages  on  the  wave  pattern,  the  DE10NA  was  also  selected  It  is  identical 
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lo  the  DEIG06  with  the  exception  that  the  appendages  are  removed.  Id  order  to  investigate  the  effect 
of  ship  length,  shorter  geosims  of  the  CVA-59  were  also  chosen:  the  CVASH  and  CVAXS  with  respec¬ 
tive  lengths  of  305  and  iw  ft  Finally,  te  order  to  investigate  the  wsversaking  of  a  strut,  which  is 
characterized  by  a  large  draft  to  beam  ratio,  SHSTNA,  the  Sharma  strut  [13]  scaled  to  a  length  of  308 
ft  was  also  selected. 

3.2  Choice  of  Exciting  Wares 

The  exciting  waves  were  taken  to  have  wavelengths  X  of  35.1,  140.4,  315.9,  and  998.3  ft.  These 
are  typical  for  Sea  States  2,  4,  6,  and  8,  respectively.  Table  2  shows  the  values  of  X,  the  corresponding 
values  of  «,  and  typical  values  of  the  amplitude  a.  The  three  largest  values  of  X  approximate  closely 
the  lengths  of  the  various  hulls,  while  the  lowest  value  of  >.  approximates  the  beam  or  draft  of  most  of 
the  hulls.  For  convenience,  the  four  particular  values  of  X  were  selected  so  that  the  corresponding 
values  of  u  coincide  precisely  with  four  of  the  30  values  for  which  force  and  motion  coefficients  are 
calculated  in  SMP. 

For  each  wavelength,  five  wave  headings  0  were  selected:  0,  45,  90,  135,  and  180  degrees.  These 
correspond  respectively  to  following,  quartering,  beam,  bow,  and  head  seas.  Thus,  a  total  of  20  dif¬ 
ferent  exciting  waves,  in  te<  ms  of  X— £  combinations,  were  used  for  each  ship. 

3.3  Computer  Program  SMPRAD 

The  above  matrix  of  computer  runs  were  made  by  using  computer  program  SMPRAD.  It  consists 
of  program  SMP  to  make  the  near  field  ^^Iculations  of  source  strength  y  (f  ,tj,{)  as  well  as  the  complex 
motion  transfer  function  Tj.  Additional  coding  was  then  developed  to  perform  the  far  field  calculations 
resulting  in  the  complex  amplitude  coefficients  A}. 


4.  FAR  FIELD  WAVE  PATTERNS 

4.1  Method  of  Presentation 

A  number  of  sample  plots  are  presented  to  give  an  indication  of  the  major  features  of  the  radia¬ 
tion  and  diffraction  wave  patterns.  The  figures  essentially  give  tht  phase  angle  <iij  (Ro,k,0,  &)  and  vari¬ 
ous  measures  of  the  magnituae  I Aj  (Ro,X,0,  /3)|  of  the  complex  amplitude  coefficient  A}  defined  in 
Eqs.  (5)  and  (6)  at  a  reference  radius  Rn  -  308  ft  This  is  equal  to  the  ship  length  for  four  of  the 
hulls.  Recall  that  j  -  1,  ....  6  refer  to  the  six  radiation  waves  due  to  ship  oscillation,  while  j  -  7 
refers  to  the  diffraction  wave  around  a  stationary  ship. 


(12) 


The  value  of  ^y  is  simply  given  by 


lk,0,  /3)  -  tan 


-i 

Im  ( Aj ) 


and  is  independent  of  R0.  The  value  of  \ A j\  is  given  by 

Uj  (k,9,  j8)i  -  -yJ(ReAj)2  +  Um  Aj)2. 


(13) 


Eqs.  (5)  and  (6)  show  that  the  vaiues  of  idyl  for  different  reference  distances  R$  and  R\,A^  and  AJX, 
are  simply  related  by 


Ujo  i  ^ 

Idyll  -J  Rq 


(14) 
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In  the  following  figures,  the  dimensional  wave  amplitude  \Aj\  is  given  in  terms  of  one  of  the  fol¬ 
lowing  dimensionless  ratios: 


Aja{k,6£) 


Idyl 


Id, 


Aj,  (k,6£)  max  [\Aj  (*,$.,£)  j,]^  (/c,fi,,0)|} 


(15a) 

(15b) 

(15c) 


where  the  subscript  a  refers  to  Aj  calculated  by  using  the  approximate  evaluation  of  Hj  given  in  Eqs. 
(8).  The  ratio  AJa  is  simply  the  amplitude  of  the  radiation  or  diffraction  wave  referenced  to  the  ampli¬ 
tude  of  the  exciting  wave.  This  is  useful  in  comparing  the  relative  wavemaking  of  the  ships  with 
respect  to  different  exciting  waves.  The  ratio  Ajt  normalizes  \Aj\  so  that  the  maximum  value  is  l. 
This  is  useful  in  comparing  the  relative  shapes  of  the  various  wave  modes  j  —  l, ..,6,7.  The  ratio  AJt  is 
similar  to  4y»  with  the  exception  that  it  is  normalized  with  respect  to  A,  using  the  accurate  and  long 
wavelength  approximations  of  H,.  This  is  useful  for  ascertaining  the  error  incurred  in  using  the  long 
wavelength  approximation  of  1^4,1. 


of 


In  all  cases,  the  wave  patterns  are  given  in  terms  of  the  variation  of  4ij  and  the  various  measures 
Aj  |  with  the  direction  6  All  of  the  curves  in  these  figures  are  drawn  with  one  of  six  line  styles 


which  are  identified  is  Curves  1  to  6.  shown  in  Figure  2. 


tJ’V'Jw-iVyi  'toJtXlXXJUSJtL  1  ftlMBCWSfati*  'JWWVWy  vfWWyiiO«> _ t  _ ■rfVVt 


p  '  4.2  rrtsMnttiro  w  Wire  P^icibs 

Figures  3a  to  3f  show  the  variation  with  9  of  the  dimensionless  amplitude  coefficient  AJa  for  the 
radiation  waves  for  the  three  vertical  and  the  three  lateral  modes.  Figures  4a  to  4f  show  corresponding 
variations  of  the  phase  angle  Due  to  the  iarge  values  of  the  lateral  mode  waves  of  the  strut-shaped 
SHSTNA  hull,  their  values  have  been  divided  by  10  in  order  tc  be  on  the  same  scale  as  the  other  hulis. 
For  the  sake  of  consistency.  Aia  for  the  vertical  modes  as  well  as  the  diffraction  waves  shown  in  Figure 
5  for  this  hull  have  also  been  divided  by  10.  This  fact  is  dearly  indicated  in  the  legend  for  these  fig¬ 
ures.  The  wavelength  X  and  heading  £  of  the  exciting  wave  considered  in  each  figure  correspond  to 
the  particular  case  out  of  the  20  x— 0  combinations  which  give  the  largest  value  of  AJa  for  the  DE1006. 

Figures  Sa  and  5b  respectively  show  the  variation  of  A1Va  and  Alta,  the  symmetrical  and  anti- 
symmetrical  components  of  the  diffraction  wave.  Figures  6a  and  6b  show  corresponding  values  for  the 
phase  angles  ifay  and  tfaL.  The  results  shown  in  these  figures  are  for  X  -  140.4  ft  and  0  *  135  deg,  a 
bow  sea.  For  four  of  the  hulls,  this  value  of  X  lies  intermediate  between  the  beam  and  length  of  the 
ship,  an  optimum  range  of  X  for  the  low  order  diffraction  theory  used  in  this  work. 

Figures  7a  to  7c  respectively  show  the  variation  of  the  normalized  amplitude  Ajb  of  the  four  verti¬ 
cal  modes  (j  -  1,3,5, 7  V)  for  X  -  315.9  ft,  0  -  135  deg,  for  the  geosims  of  different  length:  CVA-59, 
CVASH,  and  CVAXS.  Figures  8a  to  8c  show  corresponding  results  for  the  four  lateral  modes 
(J  -  2A.6.1L). 

Figures  9a  to  9d  respectively  show  ihe  variations  of  Ajb  of  the  vertical  modes  (J  -  1 .3,5,7  V)  for 
different  wavelengths  X  for  the  DEI 006.  Figures  10a  to  lOd  show  the  corresponding  variation  of  A,h 
with  X  for  the  lateral  modes  (J  -  2,4,6,71.). 

Figures  11a  and  lib  respectively  show  the  variation  of  the  diffraction  amplitude  coefficients 
A 7 ya  and  AiLa  with  wave  heading  0  for  the  previously  considered  intermediate  wavelength  x  «  140.4 
ft  (see  Figs.  5  and  6)  for  the  DE1006. 

Figures  12a  to  1 2d  respectively  show  ihe  variation  with  X  of  the  exact  and  approximate  values  of 
Aj,  of  heave,  pitch,  sway,  and  yaw  at  a  fixed  heading  of  90  deg  for  the  DE1006. 

4.3  Discussion  of  Wave  Patterns 

Leaving  aside  the  unusual  strut-shaped  hull,  SHSTNA,  Figure  3a  shows  the  expected  trend  that 
the  wave  amplitudes  for  surge  are  one  order  of  magnitude  lower  than  the  otner  radiation  modes.  Tne 
large  value  for  the  strut  case  is  that,  due  to  its  extremely  thin  waterplane  section,  there  is  little  potential 
damping  and  motions  may  be  large  near  resonance  with  resultant  large  waves.  This  also  accounts  for 
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the  larger  values  of  pitch  for  the  SHSTNA,  shown  in  Figure  3c.  Lee  [14]  discusses  thu  need  for  adding 
viscous  damping  terns  to  obtain  accurate  motion  predictions  near  resonance  for  twin  hull  ships  which 
have  strut-like  members.  Fcr  the  lateral  modes  Figures  3d  and  3e  show  that  the  wave  patterns  for  the 
SHSTNA  are  similar  to  those  of  the  other  hulls,  but  with  maximum  values  which  are  an  order  of  mag¬ 
nitude  higher.  This  is  due  to  its  large  draft.  In  the  case  of  sway,  shown  in  Figure  3f,  the  wave  ampli¬ 
tudes  for  the  SHSTNA  are  identically  zero  since  its  hull  is  symmetric  fore  and  aft  and  the  wave  heading 
is  90  deg. 

Figure  3  shows  that  the  wave  pattern  varies  significantly  with  mode  and  ship  length.  Thus,  AJa 
for  surge  and  pitch  tend  to  have  minimum  values  at  0  “  90  deg  while  the  remaining  modes  tend  to 
have  maximum  values  here,  in  both  Figures  3  and  4,  the  curves  tor  the  shortest  hull  CVAXS  tend  to 
have  the  smoothest  variation  with  9,  while  those  for  the  longest  hull  CVA-59  tend  to  have  the  sharpest 
variation.  This  is  simply  an  indication  that  a  given  value  of  X  is  longest  relative  to  the  CVAXS  and 
shortest  relative  to  the  CVA-59. 

Figures  3d  to  3f  show  that  the  removal  of  appendages  from  the  DEI 006  increases  the  maximum 
wave  amplitude  of  the  rotational  modes  of  yaw  and  roll  by  approximately  20  oercent  and  has  relatively 
little  effect  on  the  translational  sway  mode. 

In  the  case  of  the  diffraction  waves,  shown  in  Figures  5  and  6,  perhaps  the  most  striking  feature 
is  the  extremely  large  values  for  the  CVA-59,  which  has  the  longest  length.  Here,  the  exciting  wave  is 
shortest  relative  to  it,  and  hence  undergoes  the  largest  scattering.  By  the  same  reasoning,  the  CVAXS 
has  the  shortest  length  and  the  exciting  wave  undergoes  relatively  little  scattering,  resulting  in  small  dif¬ 
fraction  wave  amplitudes.  The  effect  of  hull  length  is  also  evident  in  Figures  6a  and  6b,  which  show 
that  the  phase  angle  variation  is  sharpest  for  the  CVA-59  and  smoothest  for  the  CVAXS. 

As  in  Che  case  for  the  radiation  waves,  the  behavior  of  the  strut-shaped  SHSTNA  sharply  differs 
from  that  of  the  other  hulls  of  equal  length.  It  has  the  lowest  values  for  the  vertical  diffraction  wave 
amplitude,  but  the  second  highest  (next  to  the  previously  mentioned  CVA-59  case)  lateral  diffraction 
waves.  These  are  again  due  respectively  to  its  thin  beam  and  deep  draft. 

Figures  7  and  8  show  that  the  normalized  wave  amplitudes  AJt  become  smoother  with  increasing 
ship  length  foi  the  geosim  series  CVAXS  (140  ft),  CVASh  (308  ft),  and  CVA-59  (990  ft).  The  wave 
is  fixed  at  X  -  315.9  ft  and  fZ  —  135  deg.  The  largest  change  in  the  wave  patterns  occurs  between  the 
CVAXS  and  CVASH  hulls,  and  less  of  a  change  occurs  between  the  CVASH  and  CVA-59  hulls.  Thus, 
the  shape  of  the  iong-wavelength  wave  pattern  is  generally  established  for  \=I. 


Figures  9  and  10  show  the  manner  in  which  Aj8  of  each  of  the  eight  wave  modes  approaches  the 
long  wavelength  pattern  for  increasing  wavelength  for  the  DE1006  hull  (308  ft).  The  modes  of  surge, 
pitch,  and  yaw,  which  have  minima  in  the  beam  direction  ( 9  -  90  deg),  tend  to  have  more  pronounced 
variations  with  X  than  the  remaining  modes  which  have  maxima  at  9  —  90  deg. 

Figures  11a  and  11b  show  that  the  9  -  value  of  the  maxima  of  the  dimensionless  diffraction  wave 
amplitudes  Ajy„  and  A-)La  for  the  DE1006  shifts  along  the  0  -  axis  to  always  coincide  with  the  wave 
heading  /3.  The  maximum  value  of  the  vertical  wave  A1Va  does  not  vary  with  0,  but  that  of  the  lateral 
wave  At i  varies  according  to  the  sin/3  factor  given  in  Eq.  (6).  Thus,  for  following  and  head  seas,  0  — 
0  and  180  deg.  Figure  lib  shows  the  expected  case  of  the  absence  of  lateral  diffraction  waves. 

Figures  12a  to  12d  show  the  variation  with  wavelength  X  of  the  differences  between  the  values  of 
the  normalized  amplitudes  AJt  for  four  oscillation  modes  for  the  DE1006,  using  the  exact  and  long 
wavelength  approximation  of  the  Kochin  function.  The  three  values  of  wavelength  considered  in  each 
figure  correspond  to  ratioes  k/B  —  1.0,  3.8,  and  9.1.  As  expected,  agreement  is  poor  at  k/B  —  1.0. 
However,  there  is  only  a  slight  difference  at  X/B  —  3.8,  and  nearly  identical  agreement  at  X/B  —  9.1. 
At  the  shortest  wavelength,  k/B —  1.0,  the  agreement  is  better  for  the  lateral  modes  of  sway  and  yaw 
than  for  the  vertical  modes  of  heave  and  pitch.  This  is  due  to  the  fact  that  the  long  wavelength  approx¬ 
imation  in  Eq.  (7b)  is  accurate  to  second  order  in  9  for  the  lateral  modes,  and  only  to  first  order  for  the 
vertical  modes. 

S.  SUMMARY 

The  calculation  procedure  consists  of  two  principal  parts:  the  inner  flow  and  far  field  regions. 
The  inner  flow  near  the  hull  is  calculated  by  using  the  strip-theory  Ship  Motion  Program.  The  principal 
output  quantities  of  interest  to  the  present  work  are  the  motion  transfer  coefficients  and  the  singularity 
distribution  over  the  hull  surface.  The  far  field  ship  radiation  and  diffraction  wave  patterns  are 
obtained  by  using  the  Kochin  function,  which  represents  a  weighted  integral  of  the  singularity  distribu¬ 
tion  over  the  hull  surface,  in  conjunction  with  an  asymptotic  expression  for  the  waves  due  to  an  oscil¬ 
lating  three-dimensional  source.  It  is  shown  that  a  previously  suggested  approach  of  collapsing  the  hull 
surface  singularity  distribution  to  a  line  is  equivalent  to  a  long  wavelength  approximation  of  the  Kochin 
function. 

Wave  patterns  were  calculated  for  six  ship  hulls,  ranging  in  length  from  140  to  990  ft.  Four 
wavelengths  of  the  exciting  waves  were  used,  ranging  from  35  to  998  ft.  Five  wave  headings  were  con¬ 
sidered,  ranging  from  following  to  head  seas. 
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Among  the  six  radiation  modes,  the  waves  due  to  surge  tend  to  be  an  order  of  magnitude  smaller 
than  those  of  the  other  modes.  Heave  tends  to  give  the  largest  wave  amplitudes.  The  wave  amplitudes 
due  to  surge,  pitch,  and  yaw  tend  to  have  their  minimum  values  in  the  beam  direction,  $  —  90  deg., 
while  the  remaining  three  modes  have  their  maximum  values  here.  The  major  features  of  the  long 
wavelength  wave  pattern  are  established  for  wavelengths  of  the  exciting  wave  eo.ual  to  one  ship  length. 

The  direction  of  the  maxima  of  both  the  vertical  (symmetrical)  and  lateral  (antisymmetrical)  com¬ 
ponents  of  the  diffraction  wave  coincides  with  that  of  the  wave  heading  /3.  The  amplitudes  of  the  verti¬ 
cal  component  do  not  vary  with  0  while  those  of  the  lateral  component  are  scaled  by  the  factor  sin  p. 

The  wave  amplitudes  calculated  by  using  the  long  wavelength  approximation  and  exact  evaluation 
of  the  Knchin  function  are  in  poor  agreement  for  wavelength-to-beam  ratio  \/B  of  1.0  but  is  in  excel¬ 
lent  agreement  at  >./B  ■»  3.8.  At  low  values  of  X/B,  the  approximation  is  more  accurate  for  the  lateral 
modes. 

Due  to  its  substantially  larger  draft-to-beam  ratio,  there  is  usually  a  striking  difference  between 
the  wave  patterns  of  the  strut-like  SHSTNA  hull  and  the  conventional  ship  hulls  of  comparable  length. 
The  amplitudes  of  this  hull  for  all  the  lateral  modes  are  typically  one  order  of  magnitude  larger  than 
those  of  the  other  hulis,  while  the  vertical  component  of  the  diffraction  wave  is  usually  much  smaller. 
The  case  of  the  vertical  radiation  waves  is  more  complex.  They  may  be  much  smaller  (due  to  the  nar¬ 
row  beam),  or  they  may  be  somewhat  larger  than  those  of  the  other  hulls  (due  to  large  amplitudes  of 
oscillation  near  resonance,  caused  by  small  potential  wave  damping.) 
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Table  1  —Summary  of  Hull  Geometric  Characteristics 


Ship  Hull 

L  (ft) 

B/L 

D/L 

Appendages 

1.  DEI 006 

308 

.117 

.039 

Yes 

2.  DE10NA 

308 

.117 

.039 

No 

3.  CVA59 

990 

.131 

.036 

Yes 

4.  CVASH 

308 

.131 

.036 

Yes 

5,  CVAXS 

140 

.131 

.036 

Yes 

6.  SHSTNA 

308 

.050 

.150 

No 

L  —  length,  B 

—  beam,  D  ” 

draft 

Table  2— Summary  of  Exciting  Wave  Characteristics 


Sea  State 
No. 

Wavelength 

X  (ft) 

Frequency 
«  (rad/s) 

Amplitude 
a  (ft) 

2 

35.1 

2.40 

1.0 

4 

140.4 

1.20 

3.6 

6 

315.9 

0.80 

8.2 

8 

998.3 

0.45 

22.8 

I,  -  SURGE  {,  -  HEAVE  ^  -  PiTCH 

|2  -  SWAY  $4  -  ROLL  ^  -  YAW 

Fig.  1  —  Definition  of  coordinate  system  and  six  modes  of  ship  oscillation 


Un»  Styte  1  . . *  Lrw  Sty*  3  . .  . Un»S?ytt5 


Li™  Sty*  2  -t - *•  ■  Ln.Slyt.4  ■»-  «-  Ufa  Sty*  B 


DIRECTION  (dag.) 


Fig.  2  —  Definition  of  curve  line  styles 
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DIRECTION  idag.) 

(d)  j  —  2,  swuy,  x  -  140.4  ft,  0  —  90  deg 

Fig.  3  (Continued)  —  Dimensionless  amplitude  coefficient  A,a  for  radiation  modes  for  all  six  hulls 


DIRECTION  Id-g.J 

(e1  j  »  4,  roll,  )i  -  540.4  ft.  0  -  90  de* 


oetoos 


aimrriow  id^.j 

(»)  *7^.  A  —  140.4  ft.  p  «•  125  deg 


CVAXS 


OJ  XJ  «U  BJ9  UU  HU  UB4 

DlflEETION  (d*g.) 


(b)  *H,  a  -  140  4  ft,  8  -  135  dsg 

Fig.  6  —  Phase  angles  and  <j>1L  for  diffraction  wave  for  all  six  hulls 
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DIRECTION  td*g. ) 

(b)  CVA.SH 

[.  7  ■  -  Normalized  amplitude  coefTicient  AJt  for  the  vertical  modes 
for  CVA-59,  CVASH,  and  CVAXS,  \  -  315.9  ft,  -  135  deg 
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an.a 
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aOJ  «U> 

DIRECTION 


(a)  j  “  1,  surge. 


X-34.1W 


•* - «r 


(dag.] 

0  —  13$  deg 

X»3I1W« 


-9-  XulJOMt 


X-S9BJH 


a.a 


13JJ 


1*3.3 


MU 


x.a  tcj  bj 

DIRECTION  (d*g.) 

(c)  j  —  5,  pitch.  /3  *»  135  deg 


k-36.1ft  1  *  x-315.9* 


— *•  *»i4a«e  -* — -  •••*-  x«ae.3ft 


(d)  i  ■*  71'.  diffraction,  0  ■*  135  deg 

Fig.  9  (Continued)  —  Normalized  amplitude  coefficient  Aj(,  for  the  vertical  modes  of  the  DF10O6 

for  different  waveiergths 


DIRECTION  (dag.} 

(a)  7  —  2.  sway,  p  -  90  deg 
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x.USJft 


“^k 

§ 


a 


a.0  au  hu  ianj  iso.o  ihu 

DIRECTION  (dag.) 

(b)  >  -  4,  roll,  P  —  90  deg 

Fig.  10  —  Normalized  amplitude  coefficient  AJt  for  the  lateral  modes  of  the  DE1006 

for  different  wavelengths 
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DIRECTION  tdag.l 


(c)  j  ••  6.  yaw,  $  “  45  deg 


X  =.  36.1ft  w  .  «  X=.3'5_9h 

X .  1 40.4ft  -«— - -»■  X»988J#t 


(d)  j  •*  7£,  diffraction,  0  —  135  deg 


Normalized  amplitude  coefficient  /4;J  for  the  lateral  modes  of  the  DE1006 
for  different  wavelengths 


o.o  slo  su  «ui  uo. a  mo  mo 

OIRffilUW  td*g.  j 

(a)  j  «  3,  heave 


ojj  jdo  ®o  »o  ixj:  mo  i«u 

DIS££T?«  Idsg  ! 

(b)  j  "  5,  pitch 


Fig.  12  —  Exact  and  approximate  values  of  normalized  amplitude  coefficient  Aj 
for  four  radiation  modes  of  the  DEI 006  for  different  wavelengths,  £  ”  90  deg 
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aiRECTICN  (deg.) 

(d)  j  m  6,  ytw 

Fig.  12  (Continued)  —  Exact  and  approximate  values  of  normalized  amplitude  coefficient  AJt 
for  four  radiation  modes  of  the  DE1006  for  different  wavelengths,  £  —  90  deg 
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